Enhancing the critical current in quasiperiodic pinning arrays 
below and above the matching magnetic flux 
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Quasiperiodic pinning arrays, as recently demonstrated theoretically and experimentally using a 
five-fold Penrose tiling, can lead to a significant enhancement of the critical current I c as compared 
to "traditional" regular pinning arrays. However, while regular arrays showed only a sharp peak 
in I c (&) at the matching flux $i and quasiperiodic arrays provided a much broader maximum at 

< <3>i, both types of pinning arrays turned out to be inefficient for fluxes larger than <3>i. We 
demonstrate theoretically and experimentally the enhancement of I c (&) for $ > <E>i by using non- 
Penrose quasiperiodic pinning arrays. This result is based on a qualitatively different mechanism 
of flux pinning by quasiperiodic pinning arrays and could be potentially useful for applications in 
superconducting micro-electronic devices operating in a broad range of magnetic fields. 

PACS numbers: 74.25.Wx; 61.44.Br; 74.25.Sv 



I. INTRODUCTION 

Pinning of Abrikosov vortices in superconductors is 
a topic that has attracted condsiderable interest, both 
with respect to the fundamental physical properties of 
so-called "vortex-matter" and with respect to device ap- 
plications. In the latter case, for bulk or thick film super- 
conductors, the inclusion of nanodefects, acting as pin- 
ning sites, typically randomly distributed within the su- 
perconductors, has been shown to significantly enhance 
the critical current density, which is important for in- 
creasing the current-carrying capacity of wires and tapes, 
e.g., for applications in superconducting magnets (see, 
e.g., Ref. [1[ and references therein). On the other hand, 
recent progress in the fabrication of nanostructures has 
provided a wide variety of well-controlled lithographically 
defined artificial pinning sites, and experiments can con- 
trol both the location and the strength of each pinning 
center in vortex- confining Arrays of Pinning Sites (APS) 
(see, e.g., (sMEJ ) • In this context it has been shown that 
e.g. microholes (antidots) or magnetic dots may im- 
prove the performance of micro-electronic thin film de- 
vices. Examples are the suppression of low- frequency flux 
noise in superconducting quantum interference devices 
(SQUIDs) by strategically positioned antidots 0, [?J or 
the enhancement of the quality factor of superconducting 
microwave resonators operated in magnetic fields, which 
can be achieved by incorporating antidot arrays into the 
thin film structures |8j- From a basic point of view, it 
is interesting to explore the optimum pinning landscape 
provided by APS. One convenient way to do this is to 
investigate the critical current I c versus magnetic flux <I> 
for different arrangements of pinning sites. 

Periodic artificial APS were first proposed as devices 
able of trapping magnetic flux and thus enhancing I c . 
However, the efficiency of periodic APS is strongly lim- 



ited to so-called matching fields, i.e., when the number 
of vortices is commensurate with the number of pins, 
thus resulting in strong but very narrow peaks in the 
critical current versus magnetic flux, I c (®)- This short- 
coming can be overcome by using more sophisticated 
pinning topologies, e.g., quasiperiodic (QP) APS (9l4T^. 
In particular, using a five- fold Penrose tiling (see, e.g., 
Refs. |15l-l2l|) as an APS (i.e., placing pins on the vertices 
of a Penrose tiling) results in a considerable enhancement 
of i c (^ ) ), i-e., significant broadening of the peak in I c (®), 
as was recently demonstrated theoretically [H, [T(j and 
experimentally [12|, [l3|. The underlying idea of using 
QP APS is that, contrary to regular APS, they include 
many periodicities and thus can trap vortices for mag- 
netic fluxes other than matc hing values <&i. As shown 
for Penrose-tiling APS |jj E3, EH, in addition to the 
matching- flux peak at $ = <3>i, I c (&) acquires a remark- 
ably broad maximum for <I> < 3>i. This maximum turned 
out to be even more robust with respect to variations 
of parameters of APS than the sharp peak at $ = 3>i 
HJH3]. However, in theory and experiments, I c (&) de- 
creases rapidly for fluxes <$> > 3>i for QP Penrose-tiling 
APS as well as for regular APS. Here, we discuss the es- 
sential difference in the flux pinning for magnetic fluxes 
below and above the matching value $ = <I>i and we 
demonstrate, theoretically and experimentally, how the 
critical current I c (&) can be improved in the important 
regime where $ > 3>i, by using novel artificial APS. 

The paper is organized as follows. The model and ex- 
periment are described in Sec. II and Sec. Ill, respec- 
tively. In Sec. IV, we analyze the enhancement of the 
critical current below the first matching field in differ- 
ent aperiodic pinning arrays and compare the results to 
those for a five- fold Penrose-tiling APS. The mechanisms 
of vortex "conductivity" leading to the enhancement of 
the critical current below and above the matching flux are 
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discussed in Sec. V. In Sec. V, we also present our experi- 
mental and numerical results for QP APS demonstrating 
the enhancement of the critical current for > $1. Fi- 
nally, the conclusions are given in Sec. VI. 



II. SIMULATION 

We model a three-dimensional (3D) slab, infinitely long 
in the z-direction, by a 2D (in the xy-plane) simulation 
cell with periodic boundary conditions. Note that follow- 
ing Refs. [q|, [To(| we apply periodic boundary conditions 
at the boundaries of the square simulation cell while the 
quasiperiodic array is taken smaller than the cell. This 
is done in order to prevent imposing ad hoc periodicity 
to the vortex motion in a QP APS: the free-of-pinning 
region between the QP APS and the boundary of the 
simulation cell serves as a reservoir of vortices that mim- 
ics the external applied magnetic field. This reservoir (of 
thickness of the order of few intervortex distances) erases 
the memory, i.e., the information on the coordinate of a 
vortex leaving the cell. Each additional vortex enters the 
QP sample at the current point of "least resistance to en- 
try" , similarly to the way vortices enter any other sample 
when placed in an external magnetic field. This approach 
has been successfully used in numerous simulations in the 
past (see, e.g., [1, EE [IIS ) • 

To study the dynamics of vortex motion, we per- 
form simulated annealing simulations by numerically in- 
tegrating the overdamped equations of motion (see, e.g., 
Refs. [91TIq| for a description of the method): 



cyv i rv 



e + if 
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Here ii is the total force per unit length acting on vor- 
tex z, i? v and i^ p are the forces due to vortex- vortex 
and vortex-pin interactions, respectively, if is the ther- 
mal stochastic force, and if is the driving force; is 
the velocity and rj is the viscosity. All the forces are ex- 
pressed in units of fo = <l>o/87r 2 A 3 , where $0 = hc/2e, 
and lengths (fields) are in units of A (<l>o/A 2 ). 
Following Refs. 

SEE 

we model vortex pinning by 
short-range parabolic potential wells located at positions 
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The pinning force is 
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where N p is the number of pinning sites, f p is the maxi- 
mum pinning force of each potential well, r p is the range 
of the pinning potential, 6 is the Heaviside step function, 
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The ground state of a system of moving vortices is ob- 
tained by simulating the field-cooled experiments (e.g., 
[25]). For deep short-range (5-like) potential wells, the 

energy required to depin vortices trapped by pinning sites 

(v) 

is proportional to the number of pinned vortices, Ny • 
Therefore, in this approximation, we can define the nor- 



malized critical current as follows [H, E3 : 

Ic($)=JO- 



(3) 



and study the dimensionless value J c = I c /jo, where N v 
is the total number of vortices within the simulation cell 
of area A, and jo is a constant (i.e., the current needed 
to depin N v non- interacting vortices, jo = N v jf p ). Since 
the area of the QP pattern A qp is less than A, the critical 
current Eq. (j3j) acquires a prefactor A qp /A ^0.58 foL HoL 
|26[ in the simulations below. We use narrow potential 
wells as pinning sites, with the maximum pinning force 
f p and the radius r p = 0.04 A to 0.1 A. 



III. EXPERIMENT 

Experiments were performed on d = 60 nm thick Nb 
films containing circular antidots (diameter D = 400 nm) 
at the vertices of various types of QP APS (see insets of 
Figs. 1 and 2) in an area of 200x200 /im 2 . The aver- 
age antidot density is n p ~ 0.5/im -2 , corresponding to 
a first matching field B\ = n p &o ~ 1 mT. Electric trans- 
port measurements for the determination of the critical 
current I c (&) of perforated Nb bridges were performed 
at well-controlled temperatures T and magnetic fields; 
for details of sample fabrication and experimental setup 
see Refs. [l2|, [27j. The bridges had transition tempera- 
tures T c ^ 8.4 K and transition widths varying from 15 
to 25 mK. The temperature dependence of the critical 
current at B = differed from bridge to bridge leading 
to different absolute values of I C (B = 0) at the same re- 
duced temperatures T/T c . The experimental data shown 
below, were recorded at T/T c = 0.9995, where the abso- 
lute values of I c for different samples were very similar. 
Note that according to our estimates based on a simple 
core pinning model (see Ref. [12]), the maximum value 
of the critical depinning current in quasiperiodic pinning 
arrays of antidots reached the values of about 0.5 of the 
Ginzburg-Landau depairing current, for the used param- 
eters. 



IV. CRITICAL CURRENT IN QP TILING APS 

For a five-fold Penrose-tiling APS i, E3, El El, Pi- 
ning arrays with many periodicities [28| can enhance the 
critical current [29[ for a broader range of magnetic fluxes 
as compared to periodic APS. In order to examine and 
possibly further optimize J c ($) m APS, in this section we 
study the critical current as a function of applied mag- 
netic flux, for different QP tilings. 

First, we consider a tiling consisting of the simplest 
shapes, squares and triangles (both compatible with the 
lowest-energy vortex lattices). The result for a quasiperi- 
odic square-triangle tiling [30] is shown in Fig. 1(a). For 
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FIG. 1: (Color online) Normalized critical current, J c versus normalized magnetic flux, <3>/<E>i (oc number of vortices N v ), 
for different quasiperiodic (QP) Arrays of Pinning Sites (APS) [shown by solid blue (dark grey) lines]: (a) Square-triangle 
quasiperiodic tiling, (b) Square-triangle nonpisot tiling, (c) Goodman-Strauss tiling, and (d) Nonpisot Penrose tiling. Insets 
show the corresponding tiling patterns. For comparison, the function J c ($/$i) for a five-fold Penrose-tiling APS is shown by 
a red (grey) line. The J c ($/$i) curves are shown for f P / fo — 2 and r p /X = 0.1. These values of f p and r p are also used in 
Figs. 2-4. 



comparison, Fig. 1 also shows the J c (^-dependence for 
a five- fold Penrose-tiling APS. J c (&) for the QP square- 
triangle tiling has a broad maximum at <£ m < <£>i, 
which is formed due to the re-arrangement of vortices 
between the square and triangular tiles when their den- 
sity changes. Note the absence of a peak at <3>i, indicating 
that the local matching conditions between the APS and 
the vortex lattice are not fulfilled at the given parame- 
ters of the APS (i.e., the average distance between the 
pinning sites, their maximum strength and radius). How- 



ever, a peak in J c (&) at 3>i is observed for the Penrose- 
tiling APS. Another tiling consisting of the same simple 
shapes (i.e., squares and triangles) is the nonpisot square- 
triangular tiling [30| shown in the inset of Fig. 1(b). Al- 
though these two types of tilings are characterizied by 
different inflation rules (i.e., different arrangements of 
the tiles), they contain the same elements, i.e., hexagons 
formed by six triangles surrounded by six squares. This 
similarity in the structure determines similar local match- 
ing effects between vortices and the pinning centers for 
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FIG. 2: (Color online) Flow pattern for vortices driven by an 
external current illustrating the "hole-like vortex "conductiv- 
ity" in the Shield-tiling APS for $ < $n $ « 0.95$i. The 
pinning sites are shown by filled large red (grey) circles. Vor- 
tex positions are shown for ten successive time intervals U, 
where i = 1 to 10: by open small circles, for i = 1 to 9, and 
by filled small green (grey) circles, for i = 10. The Lorentz 
force drives vortices from the left to the right. Most of the 
vortices do not move: their positions just slightly change ("po- 
larized") with respect to the centers of the pinning sites. The 
trajectories of moving vortices (from the left to the right) are 
shown by thick blue (grey) arrows. The motion of vortices re- 
sults in effective propagation of vacancies ("holes") from the 
right to the left (shown by thin dashed black arrows) , includ- 
ing "holes" entering the sample from outside, as in the region 
marked by the dotted rectangle. Newly created "holes" are 
shown by thin solid black circles, while annihilated "holes" 
(i.e., filled by vortices) are shown by thin dashed black cir- 
cles. 

these two tilings. As a result, the J c (<I>)-curve for the non- 
pisot square-triangular tiling [Fig. 1(b)] turns out to be 
practically identical to that for the quasiperiodic square- 
triangular tiling [compare to Fig. 1(a)]. An example of 
a tiling consisting of three different tiles (three types of 
rhombuses) is a Goodman-Strauss tiling [2l| shown in 
the inset of Fig. 1(c). Surprisingly, the behavior of J c (3>) 
for this APS is very similar to that of the Penrose-tiling 
APS, although it has an obvious shortcoming: it does not 
contain the maximum at the matching flux 3>i [Fig. 1(c)]. 

These examples indicate that the considered tilings 
provide, in general, a weaker enhancement of the criti- 
cal current as compared to the 5-fold Penrose-tiling APS 
which seems to be the best among these competitors. 
Let us then consider a different modification of a Pen- 
rose tiling, namely, a nonpisot Penrose tiling [2l|. This 
tiling also consists of two types of rhombuses. It is shown 
in the inset of Fig. 1(d). As shown in the main panel of 



Fig. 1(d), this APS provides an even better enhancement 
of J c ($) within nearly the same, as for the five- fold Pen- 
rose tiling, range of magnetic flux: it is characterized 
by an almost-flat plateau within the broad maximum in 
J c (3>). This improvement of the pinning properties is 
explained by the specific arrangement of the rhombuses 
in this kind of tiling, which turns out to be closer to 
the energetically favorable arrangement of vortices in the 
lattice, thus providing better vortex-pin matching condi- 
tions for various densities. For the five- fold Penrose-tiling 

aps BUS, 112,113, the maximum of the critical current 

at $ = $i is separated from the main broad maximum 
and is not as stable as the broad maximum with respect 
to changes of the parameters of the APS (i.e., the maxi- 
mum pinning strength) [s|,[iq|- The advantage of the APS 
based on the nonpisot Penrose tiling is that it provides a 
smooth transition from the matching-field configuration 
(i.e., at $ = 3>i) to lower- density pinned states (i.e., when 
some of the vertices of the small tiles become unoccupied 
[il EDI)- This provides a flattening of the J c ($)-curve 
[Fig. 1(d)] and thus increases the overall value of the 
maximum J c within the range of the plateau-like maxi- 
mum. This important improvement could be useful for 
applications of QP APS. 



V. THE "HOLE-" AND "ELECTRON-LIKE" 
VORTEX "CONDUCTIVITY" 

As can be seen in Fig. 1, all the above results obtained 
using different APS, while providing different profiles of 
J c (<!>), still have a common feature: J c drops drastically 
for $ > $i as it does in periodic APS (for both cases, 
$ > $i and $ < $i). 

Can we improve this? One might think that in practice 
we can increase the density of the pinning sites in the ar- 
ray and thus trap flux for any <I>. Indeed this might work 
to some extent, but the efficiency of an APS, especially 
QP APS, is determined not only by the ratio of the av- 
erage intervortex distance to the average period of the 
APS (<^ PS ), but also by the ratio of those distances to 
the magnetic field penetration depth A, which is specific 
for a given superconductor. To understand this, let us 
now briefly recall the essential difference between vor- 
tex lattice pinning by periodic and quasiperiodic APS. 
In the former case, if the symmetry of the vortex lattice 
coincides with that of the pinning array (e.g., for a trian- 
gular pinning array; in case of a square APS, transitions 
between the triangular, half-pinned and square lattices 
occur - see, e.g., [31]) the vortex lattice and the pinning 
array are commensurate for any <i^ v PS /A. In contrast to 
periodic APS, the pinning efficiency of a QP APS is sen- 
sitive to this ratio d^ s /X since the elastic deformations 
of the vortex lattice are involved in its collective pinning 
by QP APS. 
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FIG. 3: (Color online) Normalized critical current (experiment and theory), J c versus normalized magnetic flux, <3>/<3>i (oc 
number of vortices N v ), for different quasiperiodic (QP) Arrays of Pinning Sites (APS) [solid blue (dark grey) line]: (a) 
Wat anabe-Ito- Soma tiling, (b) Tuebingen tiling, and (c) Shield tiling. Insets show the corresponding tiling patterns. The 
function J c ($/$i) for a 5-fold Penrose-tiling APS is shown by red (grey) line, for comparison. 



This is explained by the fact that a vortex lattice is in- 
commensurate with any QP tiling (see (9l,Hoj). Thus, the 
efficiency of a QP APS does not scale with the density of 
pinning sites. On the other hand, shifting the broad max- 
imum in J c (&) towards larger values of <£, i.e., > 3>i, 
results in the loss of pinning efficiency for $ < 3>i. More- 
over, the problem of the effective pinning for magnetic 
flux for $ > $1 is of fundamental interest and, as we 
show below, is related to a qualitatively different type of 
flux pinning in aperiodic APS. 

The fast drop of J c ($) f° r $ > $i is explained by 
the enhanced mobility of interstitial vortices, which is 
rather high even in QP APS. Note that a very low con- 
centration of interstitial vortices (as compared to pinned 
vortices) is sufficient to drastically reduce the critical cur- 
rent. Using the language of conductivity in metals and 
semiconductors, we can call the motion of interstitial vor- 
tices "electron-like vortex conductivity". Correspond- 
ingly, the situation < <3>i when there are vacancies, 



i.e., unoccupied pinning sites, is similar to hole- type con- 
ductivity, since in this case the motion of vacancies is 
similar to the motion of holes in semiconductors. For 
fluxes <fr < 3>i, the concentration of "holes" is low, but 
the "vortex conductivity" is still provided by the motion 
of vortices (Fig. 2). However, in contrast to periodic ar- 
rays where this motion is possible due to the channeling 
of vortices (leading to the drop of J c (&) for < <£>i), 
the process of channeling is strongly suppressed in QP 
systems (as in metallic quasicrystals where the usual pe- 
riodic Bloch solution for electrons does not exist |32j). 
Thus this can explain the presence of the maxima and 
the fast drop of J c ($) in QP APS, for $ < $i and for 
> $i, correspondingly. 

Therefore, the problem is how to immobilize the ad- 
ditional (i.e., "electron-like") vortices for > $i. Al- 
though in a QP APS they appear to be less mobile than in 
periodic APS (no Bloch- type solutions), they can jump to 
adjacent tiles ("hopping conductivity") or "push" pinned 
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FIG. 4: (Color online) Distribution of vortices in a Shield- 
tiling APS calculated for: $ = $1 (a) and $ = 1.18$i (b). 
The vortices are shown by green (light grey) solid dots, while 
the pinning sites by red (grey) open circles. 

vortices and thus unpin them producing secondary inter- 
stitial vortices. To prevent this scenario, we propose to 
immobilize extra vortices in a special way, namely, to cre- 
ate well separated "reservoirs" for them which themselves 
are arranged in a quasiperiodic order. This will prevent 
the extra vortices from: (i) hopping-type motion (due to 
the isolation of the "reservoirs"), and (ii) channeling (due 
to the QP order). 

Among 2D tilings, the Watanabe-Ito-Soma tiling [2l|, 
[33|| (consisting of squares and rhombuses) and Tuebingen 
tiling [34] (consisting of two types of triangles) shown in 
the insets of Fig. 3(a) and (b), respectively, are possi- 
ble candidates. These tilings contain isolated "islands" 
formed by triangles and rhombuses and the density of 
pinning sites on these islands is rather high. The plots 
of J c ($) i n Figs. 3(a,b) show that indeed there is an ap- 
preciable increase of the critical current for > <l>i, al- 



though the overall value of J c for > $1 for both tilings 
is rather low. These features of J c ($) can be seen in both 
experimental and theoretical J c (&)- curves. Note that the 
experimental curves in Figs. 3(a,b,c) qualitatively (and 
to a certain extent quantitatively — see Figs. 3(a,b)) 
reproduce the features under discussion of the calcu- 
lated J c ($)- curves, i.e., the enhancement of J c ($) for 
<I> > 3>i. A more detailed quantitative comparison will 
require additional analysis of J c (3>) for different param- 
eters of pinning arrays (e.g., radii of the antidots) and 
various temperatures. Such analysis is beyond the scope 
of the present work and will be reported elsewhere [35j . 

A significant enhancement of J c (^) f° r ^ > $1 is ob- 
tained with the Shield tiling [2l|,[36| (with squares, trian- 
gles, and three-fold hexagons) [Fig. 3(c)]. The large tiles 
are excellent traps for interstitial vortices (see Fig. 4). 
These tiles are well separated (i.e., do not have com- 
mon sides) and are arranged in a QP order [see inset 
of Fig. 3(c)]. Thus, J c ($) obtains a pronounced max- 
imum of the same magnitude as the maximum at the 
first matching flux 3>i. (However the broad maximum for 

< $1 is suppressed, compared to the five- fold Penrose- 
tiling APS.) 



VI. CONCLUSIONS 

The critical depinning current J c (3>) was analyzed the- 
oretically and experimentally for different quasiperiodic- 
tiling arrays of pinning sites. We showed that a five- 
fold Penrose-tiling APS provides stronger enhancement 
of J c ($) than other 2D QP tilings for a broad range of 
the applied magnetic flux. It was demonstrated that the 
pinning properties of the Penrose-tiling APS can be im- 
proved by using a different modification of the tiling, i.e., 
nonpisot Penrose tiling. This leads to an overall higher 
value of J c for a broader range of magnetic fields. We 
proposed a new mechanism for the flux pinning in QP 
APS which implies an effective trap and space separa- 
tion of interstitial vortices. This prevents hopping-type 
motion and channeling of vortices in the APS and thus 
results in an increase of the critical current for flux values 
larger than the matching flux <3>i. These results might be 
useful for future applications in micro-electronic devices 
requiring a high J c over a broad range of the applied flux. 
Our proposal can be easily extended, mutatis mutandi, 
to other related systems, including colloidal suspensions 
interacting with pinning traps provided by arrays of op- 
tical tweezers [37] or vortices in rotating Bose-Einstein 
condensates pinned by potential landscapes created by 
co-rotating lasers [38[ . 
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